Shear and Layer Breathing Modes in Multilayer M0S2 
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We study by Raman scattering the shear and layer breathing modes in multilayer M0S2. These 
are identified by polarization measurements and symmetry analysis. Their positions change with 
the number of layers, with different scaling for odd and even layers. A chain model explains the 
results, with general applicability to any layered material, and allows one to monitor their thickness. 



The fast progress of graphene research, fuelled by the 
unique properties of this two dimensional (2d) material, 
paved the way to experiments on other 2d crystals^—. 
There are several layered materials (LMs), studied in 
the bulk since the sixties^, retaining their stability down 
to a single monolayer, and whose properties are comple- 
mentary to those of graphene. Transition metal oxides 5 - 
and metal dichalcogenides have a layered structure^. 
Atoms within each layer are held together by covalent 
bonds, while van der Waals interactions keep the lay- 
ers together 4 . LMs include a large number of systems 
with interesting properties^. E.g., NiTe2 and VSe2 are 
semi-metals 4 , WS 2 6 , WSe 2 7 , MoS 2 8 , MoSe 2 , MoTe 2 , 
TaS 2 9 , RhTe 2 , PdTe 2 are semiconductors 4 , h-BN, and 
HfS 2 are insulators, NbS 2 , NbSe^, NbTe 2 , and TaSe 2 
are superconductors^, Bi 2 Se3«ii, Bi 2 Te3ri show thermo- 
electric properties^ and may be topological insulators 1 ^. 
Similar to graphite and graphene, the LM properties 
are a function of the number of layers (N). The com- 
binations of such 2d crystals in 3d stacks could offer 
huge opportunities in designing the functionalities of such 
heterostructures 1 -! 2 -. One could combine conductive, in- 
sulating, superconducting and magnetic 2d materials in 
one stack with atomic precision, fine-tuning the perfor- 
mance of the resulting material 1 , the functionality being 
embedded in the design of such heterostructuresi 

Amongst these LMs, MoS 2 is a subject of in- 
tense research because of its electronic* 1 ^ and op- 
tical properties 1 ^, such as strong photoluminescence 
(PL) 14,15 , electroluminescence 1 ^, controllable valley and 
spin polarization 1 ^"—. A single layer MoS 2 (1L- 
MoS 2 ) consists of two planes of hexagonally arranged 
S atoms linked to a hexagonal plane of Mo atoms 
via covalent bond s 14 ! 20 " — . In the bulk, individual 
MoS 2 layers are held together by weak van der Waals 
forces 2 ^ - — . This property has been exploited in lubri- 
cation technology^ and, more recently, enabled the iso- 
lation of lL-MoS^i 5 -^. While bulk MoS 2 is a semi- 
conductor with a 1.3eV indirect band gap2£, 1L-MoS 2 
has a 1.8eV direct band gap 14 i 15 . The absence of in- 
terlayer coupling of electronic states at the T point of 
the Brillouin zone in 1L-MoS 2 15 i 27 results in strong ab- 
sorption and PL bands at-1.8eV (680nm)iii£. 1L-MoS 2 
field effect transistors (FETs) show both unipolar— and 
ambipolar— transport, with mobilities>500cm 2 F _1 s _1 



and on-off ratios up to 10 9 [29ll30| . !L-MoS 2 is also a 



promising candidate for novel optoelectronic devices 1 ^, 
such as photo detectors 3 -^— and light-emitting devices 
operating in the visible range. 

Raman spectroscopy is the prime non-destructive char- 
acterization tool for carbon material s 34 ! 35 , in particular 
graphite^ 5 -—, singl o 35 ! 39 and multilaye r 35 ! 39 graphene. 
The Raman spectrum of graphene consists of two fun- 
damentally different sets of peaks. Those, such as D, 
G, 2D, etc, due to in-plane vibrations 3 ^—, and others, 
such as the shear (C) modes 40 and the layer breath- 
ing (LB) modes (LBMs) 41 43 , due to relative motions 
of the planes themselves, either perpendicular or par- 
allel to their normal. Albeit being an in-plane mode, 
the 2D peak is sensitive to N since the resonant Ra- 
man mechanism that gives rise to it is closely linked 
to the details of the electronic band structur o 35 i 39 i 44 , 
the latter changing with N 45 i 46 , and the layers relative 
orientation 4 ^. On the other hand, the C modes and 
LBMs are a direct probe of Ni2r— , since the vibra- 
tions themselves are out of plane, thus directly sensitive 
to N. The success of Raman scattering in characteriz- 
ing graphene prompted the community to extend this 
technique to other LMs, from bulk to monolaye r 40 ! 47 " — . 
E.g., the Raman spectrum of bulk MoS 2 consists of two 
main peaks~382, 407cm _1 153.51]. These are assigned 
to E\ g (in-plane vibration) and A\ g (out of plane vibra- 
tion) mode a 53 i 54 . The E\ g mode red-shifts, while the 
A\ g mode blue shifts with increasing N 49 i 55 . The E\ g 
and Ai g modes have opposite trends when going from 
bulk MoS 2 to 1L-MoS 2 , so that their difference can be 
used to monitor N 4 ^. The E\ g shift with N may be at- 
tributed to stacking-induced structure changes or long- 
range Coulombic interlayer interaction s 49 ! 55 , while the 
Ai^ shift is due to increasing restoring forces as addi- 
tional layers are adde d 49 ! 55 , however, further work is still 
needed to fully clarify and assign these trends, but this 
is not the subject of the present investigation. 

Instead, our focus here is on the C and LB modes 
that appear in the low frequency region in the various 
LMs 56 . These have been extensively studied in multilayer 
graphene^ - 42 1 57 . Unlike graphite and graphene, most 
LMs consist of more than one atomic element. E.g., each 
MoS 2 layer contains one Mo plane sandwiched by two S 
planes, while Bi 2 Se3 contains two Bi and three Se planes. 
This makes their lattice dynamics more complex than 
multilayer graphene, starting from the symmetry and 
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FIG. 1: (a) Lattice structure and linear chain models of 2 and 3L-M0S2. Inversion symmetry applies to 2L, not 3L-M0S2. (b) 
Shear modes and (c) LBMs in 2 and 3L-M0S2. (d) High frequency optical vibration modes for 1 and 2L-M0S2. The symbol 
under each mode is its irreducible representation. R or IR indicate if the mode is Raman, or Infrared active, or both. 



force constants. Even NL-M0S2 belong to point group 
Dsh with inversion symmetry, while odd NL-M0S2 corre- 
spond to D^h without inversion symmetry 5 -^. There are a 
few reports on C and LBMs in LMs other than graphene. 
Ref. [59( reported them in bulk samples as:— 21.5cm- 1 
(C),32.5cm- 1 (C), -SOcm-^LBM) for As 2 Se 3 at 15K; 
-27cm- 1 (C),38cm- 1 (C), ~60cm _1 (LBM) for As 2 S 3 
at 15K; -34cm- 1 (C), -56cm- 1 (LBM) for MoS 2 ; 
-22cm- 1 (C) for GaS; -56cm- 1 (C) for GaSe. Only 
Refs (5lll52j reported some of these for non-bulk samples. 
In particular, Ref.[5l| studied the Raman spectrum of 
one shear mode for 2, 3, 5, 6 and IOL-M0S2. Ref.[52| ob- 
served only one set of C and LB modes for 2 to 6L-M0S2 
and 9L-M0S2. Both did not consider the symmetry differ- 
ence between odd- and even-NL M0S2, e.g. they assigned 
the C mode in odd-NL M0S2 as E\ , but, as we show 

later, this is instead E' . Ref.(5lj suggested that the scal- 
ing rule of the C mode in multilayer graphene^ cannot 
be extended to few layer M0S2, opposite to the results 
presented in Ref. [52]. Ref. |52| wrote that LBMs scale 
as 1/N, as predicted by Ref.|60j with an assumption of 
strong coupling between layers and substrate. However, 
it is not clear whether such strong coupling actually ex- 
ists. Furthermore, even though LBMs are optical modes, 
an acoustic atomic displacement for such modes was pre- 
sented in Fig. lb of Ref. (52| , with no symmetry analysis. 
Therefore, all symmetries, force constants, possible role 
of interactions between layers and substrate, and mode 
scaling with N still need to be fully understood. More 
importantly, it would be desirable to establish a general 
model to describe the evolution of C and LB modes with 
N in any LMs, not just M0S2. 

Here we measure the shear and layer breathing modes 
for NL-M0S2 up to 19L-MoS 2 , and bulk MoS 2 . We iden- 
tify several groups of modes with frequencies dependent 
on N. Samples with even and odd N show different scaling 
laws with N, due to their different symmetry. A simple 



chain model can account for the observed trends, and can 
be extended to other LMs. 

NL-M0S2 samples are produced from bulk M0S2 (SPI 
Supplies) by mechanical exfoliation, following a similar 
procedure to that used for graphene layers^Si. NL- 
M0S2 are supported on a Si wafer with 93nm Si02, 
which is used as substrate in order to make the sam- 
ples visible by eye. The layer thickness is determined by 
optical contrast 62 and atomic force microscopy 55 . Ra- 
man measurements are performed using a Jobin-Yvon 
HR800 system equipped with a liquid nitrogen cooled 
charge-coupled detector. The excitation wavelength 
is 532nm from a diode-pumped solid-state laser. A 
power— 0.23mW is used to avoid sample heating. The 
laser plasma lines are removed using a BragGrate band- 
pass filter (OptiGrate Corp), as these would appear 
in the same spectral range as the modes of interest. 
The Rayleigh line is suppressed using four BragGrate 
notch filters with an optical density 3 and a spectral 
bandwidth— 5- lOcm -1 . This configuration is similar to 
that used in Ref.[Io| for multi-layer graphene. The spec- 
tral resolution is— 0.6cm -1 , as estimated from the the 
Rayleigh peak full- width at half-maximum (FWHM). 

Bulk M0S2 and 2L-M0S2 belong to the space group 
P6s/ mmc (point group De/i) 6 - 3 , with unit cell consisting 
of two Mo atoms in sites with point group Da^, and four 
S atoms in sites with point group C3^, as shown in 
FigJIJi. There are 18 normal vibration modes 63 . The 
factor group of bulk and 2L-M0S2 at T is De^, the same 
as the point group^i. The atoms site groups are a sub- 
group of the crystal factor group^. The correlation 5 - 3 - of 
the Mo site group D^h^ S site group Cs v ^ and factor group 
D^h allows one to derive the following irreducible repre- 
sentations for the 18 normal vibration modes at r 53 i 65 : 
T= Ai g + 2A 2u + 2B 2g + Bi u + E\ g + 2Ei u + 2E 2g + E 2u , 
where A 2u and E\ u are translational acoustic modes, 
Ai g , E\ g and E 2g are Raman active, A 2u and E\ u are 
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FIG. 2: (a) Raman spectra of 5, 6L and bulk M0S2. (b) Raman spectra of 5 and 6L-M0S2 measured for XX (red) and YX 
(blue) polarizations. The irreducible representation of each mode is indicated. 



infrared (IR) active. The E\ g and A\ g modes and one 
of the doubly degenerate E23 modes, E^, as shown in 
FigJTJi for 2L and bulk M0S2, give rise to Raman modes 
above 200cm- 1 [H|. Only A lg (-408cm- 1 in bulk and 
-405cm- 1 in 2L-MoS 2 ) and (-382cm- 1 in bulk and 
—383cm- 1 in 2L-M0S2) can be observed when the laser 
excitation is normal to the sample basal plane 49 . The 
other doubly degenerate E 2g mode, E^, and one B 2g 
mode, B^, are shear and LB mode a 53 i 54 i 66 . E 2g cor- 
responds to a rigid-layer displacement perpendicular to 
the c axis (C modes), while B 2g corresponds to rigid-layer 
displacements parallel to the c axis (LBMs), as shown in 
Fig©D,c for 2L and bulk MoS 2 . 

IL-M0S2 has Dsh symmetry, with three atoms per unit 
cell££. The irreducible representation of Ds h 63]6b gives: 
r= 2A 2 J rA l J r 2E +E , with A2 and E acoustic modes, 
A 2 IR active, A x and E Raman active, and the other 
E both Raman and IR active. The E and A x modes, 
FigJTJi, were previously detected in the Raman spectra of 
IL-M0S2 at-384 and-403cm- 1 [i9|. Of course, no rigid- 
layer vibrations can exist in IL-M0S2. 

3L-M0S2 has the same point group (Dsh) as IL-M0S2, 
with A 1 and A 2 corresponding to LBMs, and E and E 
being shear modes, Fig^T]3,c. Systems with even N belong 
to point group (with inversion symmetry), while odd 
N correspond to D^h (without inversion symmetry)^. 



Therefore, it is convenient to denote each mode of NL- 
M0S2 by the corresponding irreducible representation ac- 
cording to the their point group, and then determine if 
they are Raman, IR active, or inactive. 

NL-M0S2 has 9N-3 optical modes: 3N-1 are vibrations 
along the c axis, and 3N-1 are doubly degenerate in- 
plane vibrations. For rigid-layer vibrations, there are N-l 
LBMs along the c axis, and N-l doubly degenerate shear 
modes perpendicular to it. When N is even, there are 
Raman active LBMs and ^ doubly degenerate shear 
modes. When N is odd, there are LBMs and N-l 
doubly degenerate shear modes. The inter-layer distance 
in LMs is much larger than the in-plane bond length, 
e.g. in M0S2 the inter-layer distance is— 6.7A, while the 
in-plane bond length is— 3.2A& Thus, the in-plane opti- 
cal modes may not strongly depend on N. However, the 
interlayer coupling dominates the lattice dynamics of the 
rigid-layer vibrations, so that LB and shear modes will be 
very sensitive to N. E.g., Figj2^i shows the Raman spectra 
of 5L, 6L and bulk M0S2. In the high frequency region 
above 200cm- 1 , the E\ g (-384cm- 1 ), A lg (-409cm- 1 ), 
2LA(M) (-453cm- 1 ) and A 2u (-463cm- 1 ) modes are 
detected both in bulk and 6L-M0S2. Although the no- 
tation in 5L-M0S2 is different from 6L and bulk M0S2 
because of the crystal symmetry, the modes (E , A x ) are 
also observed in 5L-M0S2. The lineshape and peak posi- 
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FIG. 3: (a) Stokes and anti-Stokes Raman spectra of ONL-M0S2 in the low frequency range, (b) Stokes and anti-Stokes Raman 
spectra of ENL-M0S2. The spectrum of bulk M0S2 is also included in (a,b). Dashed and dotted lines in (a,b) are guides to the 
eye. (c) Position of typical C and LB modes as a function of N. (d) FWHM of C and LBM as a function of N. Solid lines in 
(c,d) are guides to the eye 



tions in the high frequency region for 5L and 6L-M0S2 is 
almost identical, both being similar to bulk M0S2. In the 
low frequency region below 100cm -1 , there is only one 
Raman peak^33cm -1 , i.e., E\ g , in bulk M0S2. However, 
as discussed above, there should exist 6 Raman active 
modes for 5L and 3 for 6L-M0S2. Of these, 4 and 3 shear 
modes should be doubly degenerate for 5L and 6L-M0S2, 
respectively. Experimentally, we observe 3 modes below 
60 cm -1 in 5L and 4 in 6L-M0S2, as shown in Figj2ji. 

The LB (A ± ) and shear (E , E ) Raman tensors in odd 
NL-M0S2 (ONL-M0S2) and shear (E$ g ) Raman tensor in 
even NL MoS 2 (ENL-M0S2) are" 
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E (shear, ONL) : 



E (shear, ONL) : 



El g (shear,ENL) : 



We do not discuss the LBM (B% g ) in ENL-MoS 2 since it is 
Raman inactive. These tensors show that, in backscatter- 
ing, the A 1 modes in 5L-M0S2 should appear only under 



unpolarized XX configuration, and E should exist under 
both unpolarized XX and polarized YX, while E should 
not appear for either XX or YX. Here XY indicates two 
mutually perpendicular axes within the basal plane of 
NL-M0S2, the first being the polarization direction of 
the incident laser, the second the analyzer's polarization. 
For 6L-M0S2 under back-scattering, the E\ g modes exist 
for both XX and YX configurations. 

In FigJSJ} two sharp peaks are observed under both 
XX and YX configurations at~19 and^30cm -1 for 5L, 
and~23 and~32cm -1 for 6L-M0S2. According to the 
symmetry analysis discussed above, we assign these to 
E in 5L and E\ g in 6L-M0S2. Two broad peaks are 
observed for XX measurements at~17 and 47cm -1 for 
5L-M0S2, which we assign to A 1 . Note that the lower E 
mode of 5L-M0S2 at~19cm -1 cannot be fully resolved for 
XX measurements due to the presence of the broad A 1 
mode~17cm -1 . We also detect two Raman modes~15 
and 41cm -1 in 6L-M0S2 consistent with what should be 
optically silent B\ LBMs, as discussed later. The sym- 
metry, polarization and position are summarized in Ta- 
ble I for the shear and LB modes of 5 and 6L-M0S2. 
Although the in-plane modes in 5 and 6L-M0S2 above 
200cm -1 are almost identical in frequency and lineshape, 
the shear and LB mode positions below 100cm -1 are dif- 
ferent. The frequencies of all 5L-M0S2 LBMs are higher 
than in 6L-M0S2, while all shear modes are lower. 

FigSi,b show the low frequency Raman measurements 
for NL-M0S2, with N=1...19, as well as bulk MoS 2 . Since 
the point group of ONL-M0S2 (Dsn) is different from 
ENL-M0S2 (Dqh)i we plot the Raman spectra of ONL- 
(FigEk) and ENL- (FigEb) MoS 2 in two panels. Bulk 
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FIG. 4: Position of shear modes as a function of N for (a) ENL- and (b) ONL-MoS2- Position of LBMs as a function of N for 
(c) ENL- and (d) ONL-M0S2. The blue open circles are the experiment data. The red crosses are the diatomic chain model 
calculations. The irreducible representation of each mode is also indicated. Solid lines are guides to the eye. 



TABLE I: Symmetry, polarization and experimental (exp.) 
positions of shear and LB modes in 5 and 6L-M0S2 
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MoS 2 is included both in ONL-M0S2 and ENL-M0S2 
panels because we cannot distinguish its parity in a (non- 
infinite) bulk sample. Of course there are no LB nor 
shear modes in IL-M0S2, as confirmed in Figj3^i. The two 
spikes~4.55cm -1 , with weaker intensity for thicker M0S2 
flakes, are due to Brillouin scattering of the LA mode 
from the Si substrate^. This is confirmed by determin- 



ing the elastic constant en from ch=p^ 2 7t 2 /(^ 2 +^ 2 )/cq, 
where p is the Si density, v is the LA mode frequency, 
ko=27r/Ao, Ao is the incident light wavelength and (77+ift) 
is the Si complex refractive index. The en determined 
from our Raman measurements (~1.65x 10 n Pa) is con- 
sistent with 1.66xlO n Pa measured by ultrasonic wave 
propagation^!. In 2 and 3L-M0S2, two Raman peaks are 
observed, with two peaks overlapping in 3L-M0S2. More 
Raman peaks are observed for thicker M0S2 flakes. 

We classify all low frequency Raman peaks into two 
categories. Those that stiffen for increasing N, linked 
by dashed lines, and those softening with N, linked by 
dotted lines. One shear mode stiffens from~22.6cm -1 
in 2L to 32.5cm -1 in 19L-MoS 2 , while one LBM softens 
from-40.1cm -1 in 2L to 4.7cm -1 in 19L-MoS 2 , FigEt. 
The two sets of modes have different FWHM, Figj3]i. In 
2L-M0S2, FWHM(LBM-40.1cm -1 ) is-9.6cm -1 , much 
larger than that (~0.8cm -1 ) of the shear mode at 
22.6cm -1 . All data for ENL- and ONL-M0S2 are sum- 
marized in Fig^i-d. According to group analysis, there 
should be no Raman active LBMs in ENL-M0S2. How- 
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ever, a set of peaks are observed in ENL-M0S2, with 
the same polarization behavior as LBMs in ONL-M0S2, 
e.g. the two Raman modes~15 and 41cm -1 in 6L-M0S2 
in FigJJJ). Because their measured positions match well 
those predicted for LBMs, i.e. the B\ g modes in ENLs 
as discussed later, they are included in FigHh. 

Since a M0S2 layer consists of two types of atoms, S 
and Mo, we implement a diatomic chain model (DCM) to 
explain the data. FigJIJi shows the ball and stick model 
for 2 and 3L-M0S2. Only two force constants are needed 
to describe the vibrations: a ss and a sm , where a ss is 
the force constant per unit area between two nearest S 
planes in two adjacent layers, and a sm that between the 
nearest S and Mo planes within a M0S2 layer. Their com- 
ponents perpendicular to the basal plane, ajr s and ajr m1 
determine the LBMs lattice dynamics, while those paral- 
lel to the basal plane, a\ s and aim-, determine the shear 
modes dynamics. The reduced mass for a S (Mo) plane, 
ms (^Mo), is its atomic mass per unit area. In M0S2, 
rris = 0.6 x 10 -7 g/cm 2 and rriMo = 1-8 x 10 -7 g/cm 2 . 
ajr m and aim can be estimated from the high frequency 
A 1 and E modes of IL-M0S2, for which 9x9 dynami- 
cal matrices can be const ructed a nd solved analytically. 
We get uo A ' = (l/27rc)>/2a^ n / /i, with \i = 2ms- The 
atom displacement eigenvectors show that the vibration 
directions of the two external S atoms are opposite along 
the c axis, while the center Mo atom stays still, as shown 
in Fig lUi for IL-M0S2, corresponding to a spring con- 
nected by two S atoms with a force constant per unit area 
2a^- m , since the Mo atom stays still at the spring equilib- 
rium position. We measure uo A ' ~403cm -1 in IL-M0S2. 

This gives ajr m = 3.46 x 10 21 7V/m 3 . We also get uj e > = 

(l/27rc)^2aL//i, where l//i = l/m Mo + l/(2m 5 ). The 
atom displacement eigenvectors indicate that the vibra- 
tion directions of the two S atoms are opposite to the cen- 
ter Mo atom, along the basal plane, as shown in FigJT]i 
for IL-M0S2, corresponding to a spring connected by two 
S atoms and one Mo atom with a force constant per 
unit area 2aJ m . From the experimental 384cm -1 , we get 
aim = 1.88 x 10 21 N/m 3 . Note that the weak interaction 
of the two S planes in the M0S2 layers is not included 
because the S-S plane distance is twice the S-Mo one. 

To understand the shear modes of NL-M0S2, the layer 
coupling between two nearest S planes in the two adjacent 
layers should be included. 3Nx3N dynamical matrices 
can be constructed for NL-M0S2. By numerically solv- 
ing the eigen-equation for NL-M0S2, we get the eigenfre- 
quencies and corresponding eigenvectors. By fitting these 
to our experimental data we get: 0^=8. 90 x 10 19 7V/m 3 

and aL=2.82 x 10 19 7V/m 3 . Multiplying aj- s and al s 
by the unit cell area gives the interlayer force con- 
stants, ^ ear =2.5N/m and /^/ M =7.8N/m. They agree 
well with those for bulk samples reported in Ref.|59j 
(kf hear =2. 7N/m and k c 1 omp =7A'N/m), derived by consid- 
ering S-Mo-S as a rigid-laye r mass unit, and deducing the 
force constants from uo = y/ki/fii, with uj the rigid-layer 



frequency and [i\ the reduced mass. Multiplying a\ s by 
the equilibrium distance between two adjacent M0S2 lay- 
ers gives a shear modulus~18.9GPa, in good agreement 
with that measured for bulk MoS^, from phonon disper- 
sion curves determined by neutron scattering, and X-ray 
measurements of the linear compressibilities. 

The eigenvectors of the rigid-layer vibrations in 2 and 
3L-M0S2 derived from the corresponding eigen-equations 
are depicted in Figd3,c. Applying symmetry analysis 
to the corresponding NL-M0S2 eigenvectors, we assign 
the irreducible representations of the corresponding point 
group to each mode. The eigenfrequencies of Raman 
active rigid-layer vibrations for E\ g (C modes in ENL- 

MoS 2 ), E (C modes in ONL-M0S2), and Ai (LBMs 
in ONL-M0S2) are summarized in FigsSk,b,d, respec- 
tively. The eigenfrequencies of the Raman inactive B\ 
(LBMs) in ENL-M0S2 are also included in Figgt. As 
illustrated in FigHJ the model calculations are in good 
agreement with experiments, including the Raman inac- 
tive B\ g . This suggests that the Raman inactive LBMs 
(B^g) in ENL-M0S2 might be observed, with polarization 
behavior identical to the A\ (LBMs) in ONL-MoS 2 . 

We now consider the evolution of the rigid-layer vi- 
brations with increasing N based on symmetry analysis. 
In FigS^c, one C mode (E^g) and one LBM (B^g) are 
observed in 2L-M0S2. Each splits in two branches with 
increasing N, one stiffening, the other softening with N. 
A new mode appears when N increases up to 4N + 2, 
TV = 1,2,3,..., and it splits into two branches again for 
higher N. The C modes (E) (Figgb) and LBMs (A[) 
(FigHJi) in ONL-M0S2 exhibit similar trends with N as 
the C modes (^) (Fig Hi) and LBMs (£^)(Figgfc) in 

ENL-M0S2, but with decreasing frequency for E and in- 
creasing for A 1 . Connecting each branch with solid lines 
shows that these form series of cone-like curves, FigH^i-d. 
The number of LB and shear modes in ONL- and ENL- 
M0S2 increases with N. However, in the experiment, no 
more than 3 of them are observed. In both ONL- and 
ENL-M0S2, most of the observed shear modes are from 
the upper branch, and their frequencies stiffen with in- 
creasing N, while most of the LBMs are from the lower 
branch, and their frequencies soften with increasing N. 

FigE] plots the positions of all the observed and 
calculated rigid-layer vibration modes. The shear 
mode at 22.6cm -1 (E% g ) in 2L-MoS 2 blue-shifts to 
28cm -1 (£; / ) in 3L-MoS 2 , and to 32.5cm -1 (£; / ) in 19L- 
MoS 2 , reaching~32.7cm -1 (E% g ) in bulk MoS 2 . In multi- 
layer graphene 40 , the ratio of C peak positions in bulk 
and 2LG is uJbuik/^2LG = V2. In our M0S2 measure- 
ments we have 32.7/22.6 = 1.447, very close to y/2. On 
the other hand, the LBM at 40cm -1 (B% g ) in 2L red- 
shifts to 29cm -1 (A / 1 ) in 3L-MoS 2 , and Scm -1 ^) in 
I9L-M0S2. The blue-shifted branch reaches 56.8cm -1 
in 19L-MoS 2 , close to the LBM value in bulk MoS 2 73 . 
The relative displacements in 2L-M0S2 between Mo and 
two S atoms within one plane for rigid-layer vibra- 
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FIG. 5: Position of (a) shear modes and (b) LBMs as a func- 
tion N. The blue open circles are the experimental data. The 
diameter of the circles represents the Raman intensity of each 
mode. The red crosses are calculations based on the diatomic 
chain model. The black solid lines and gray dashed lines 
are, respectively, fitted by cj(N) = cj(2No)y/l + cos(Nott/N) 
(N>2N ) and u(N) = u(2N )^/l - cos(N 7r/N) (N>2N ) 
for the branches originating from 2L-(No=l), 6L-(No=3) and 
10L(No=5)-MoS2 based on the monatomic chain model 

tion modes are around 0.6%, decreasing with increas- 
ing N. We can thus further simplify the model col- 
lapsing an entire layer in a single ball. So we con- 
sider a reduced monatomic chain model (MCM). Tak- 
ing one layer as a ball with mass (rriMo + 2ra s ) and 
interlayer bonding ajr s for LBMs, and al s for shear 
modes, we get: uj L bm (l/V^c) \f ai s /(m M o + 2m s ), 

uj c = (1/ V^ttc) yj als/ (ttimo + 2m s ) for 2L-M0S2, with 
^lbm the LBM position, and uoc the C peak posi- 
tion. The corresponding uolbm (40.8cm -1 ) and uoc 
(23.0cm _1 ) are in good agreement with those from the 
DCM (uj lbm =40.3 cm -1 and cj c =22.9cm -1 ) and the ex- 
perimental data (ojlbm =/ ^0.2 cm -1 and cjc=22.6cm -1 ). 
We can now solve the eigen-equation analytically and 
find the relation between position and N both for shear 
and LB modes. These vibration modes can be assigned 
to several branches, as shown in FigEl A new branch 
will emerge from each ENL-M0S2, i.e. 2,4,6,8L...-MoS 2 , 
at about the same position as that of the C mode or 



LBM in 2L-M0S2, then splitting into two subbranches, 
one blue-shifting, the other red-shifting with increas- 
ing N. The observed C modes are usually in the high- 
frequency subbranch, while the corresponding LBMs are 
usually in the low- frequency one; these are connected 
by solid lines in Figj5] For the branches originating 
from each ENL-M0S2, the frequency as a function of 
N is u(N) = cj(2iVo) v 7 ! + cos(N tt/N) (N>2N , with 
N an integer: 1,2,3,4....) for the high-frequency sub- 
branch, and u(N) = o;(2AT ) A /i - cos(N Q 7r/N) (N>2N ) 
for the low- frequency one. E.g., for the high-frequency 
C subbranch originating from 2L-M0S2, we have No=l, 
thus u c (N) = uj c (2)yjl + cos(tt/N) (N>2), where 
co>c(2)=23.0 cm -1 . If we replace ooc(2) ~23.0cm -1 in 
M0S2 with cjc(2) ~31cm -1 , this relation describes the 
C peaks in multilayer graphenes^. Similarly, ujc(N) = 
ooc(6)y/l + cos(3tt/N) (N>6) for the high-frequency 
subbranch originating from 6L-M0S2, and uoc(N) = 
^c(10)\/l + cos(57r/iV) (N>10) for the high-frequency 
subbranch originating from IOL-M0S2, where co>c(6) and 
cjc(10) are almost the same as uoc(2). The observed 
LBMs in M0S2 mainly come from the low- frequency sub- 
branches. The relation between frequency and N in the 
low- frequency subbranch originati ng from 2, 6 an d 10L- 
MoS 2 are: uo LBM (N) = u) L BM (2)y/l - cos (tt/N) (N>2), 
ulbm(N) = u LB M(6) \/l-cos(37r/N ) (N>6) and 
ulbm(N) = ou L bm(10)^1-cos(57t/N) (N>10), where 
^L J BM(2)=40.8cm -1 . Also, ujlbm(6) and ujlbm(10) are 
almost the same as ojlbm(2). These relations match well 
the experiments, Fig^5^,b. Fig j5] also shows another sub- 
branch for 2, 6 and IOL-M0S2 both for shear and LB 
modes indicated by gray dashed lines. Only one or two 
modes are detected for these subbranches, with positions 
in good agreement with the model predictions. 

Note that any coupling between supported M0S2 and 
the substrate is not included in our chain models. The 
excellent agreement between experiments and model pre- 
dictions means that the coupling between M0S2 and the 
substrate does not play a major role, the scaling with 
N being only determined by the interaction between 
the M0S2 layers. Indeed, for the suspended multilayer 
graphene in Ref.[io|, no coupling was considered, and 
the C scaling with N was also well described by a MCM. 

In principle, our chain model can be extended to 
predict rigid-layer vibrations in other LMs. The gen- 
eral approach is to calculate the reduced mass for the 
monolayer of a given material, and then measure C and 
LBMs in 2L samples. One can then predict the rela- 
tion between frequency and N for the different branches 
in any LM. E.g., the theoretical positions of the C 
and LB modes in 2L-hBN are~38.6cm -1 and 85.6cm -1 , 
respectively^. Our model predicts that the C mo de gen- 
erates two branches, ou(N) = 38.6^1 + cos (tt/N) (N>2) 
at higher frequency, and ou(N) = 38.6^1 — cos(tt/N) 
(N>2) at lower. The LBM also generates two branches, 
cu(N) = 85.6^1 + cos(ir/N) (N>2) at higher-frequency, 
and uj(N) = 85.6^1 - cos(tt/N) (N>2) at lower. Simi- 
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larly, the C mode (-38.6cm- 1 ) and LBM (-85.6cm- 1 ) 
in 4L-hBN, will also generate two branches, and so on. 
Thus, we can predict all the C and LB modes in NL-hBN. 

In conclusion, we characterized single and few layer 
M0S2 by Raman spectroscopy. We observed rigid-layer 
vibrations both for shear and layer breathing modes, 
assigned to the irreducible representations of the point 
group which the sample belongs to, as confirmed by po- 
larized Raman spectroscopy. These change with number 
of layers, with different scaling for odd and even layers. 
A diatomic chain model, combined with group theory, 



explains the observed trends. Furthermore, a reduced 
monatomic chain model can be used to describe the shear 
and layer breathing modes in M0S2 and any other layered 
material with any number of layers. 
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